Quantum chaotic patterns in the E®(&i + 62) Jahn- Teller model 
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We study statistical properties of excited levels of the E®(J>i +62) Jahn- Teller model. The 
multitude of avoided crossings of energy levels is generally claimed to be a testimony of quantum 
chaos. We found that apart from two limiting cases (E CS> e and Holstein model) the distribution 
of nearest-neighbor spacings is rather stable as to the change of parameters and different from the 
Wigner one. This limiting distribution assumably shows scaling ~ y/S at small S and resembles 
the semi-Poisson law P(S) = 4S'exp(— 25) at S > 1. The latter is believed to be universal and 
characteristic, e.g., at the transition between metal and insulator phases. 

PACS numbers: 05.45.-a,31.30.-i,63.22.+m 



Phonon spectra of two-level electron systems - two 
phonon Exe Jahn- Teller (JT) model with rotation sym- 
metry and one-phonon exciton models HHSHIE0 
show up remarkable features: multiple avoided level 
crossings (MAC) and localized ("exotic") excited states 
at certain quantum numbers. These phenomena are typ- 
ical for chaotic spectra (e.g., in complex nuclei) usually 
associated with underlying non-integrable Hamiltonians 
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The reflection symmetric two-level Hamiltonians con- 
tain a hidden nonlinearity due to the phonon assis- 
tance of the tunneling term which reveals explicitly af- 
ter the elimination of the electron degrees of freedom 
El E3- 0|- Semiclassical approaches to the 
phonon dynamics handle this nonlinearity in different 
ways. One or another decoupling method results in los- 
ing different amounts of the quantum information, and 
hence to controversial conclusions about possible classical 
chaotic behavior Q. 

The adiabatic approach to the quantum Exe JT 
model [rotation symmetric version with two vibron (bo- 
son) modes, one symmetric and the other antisymmet- 
ric against the reflection] applies in a limited range of 
validity for the strong electron-phonon coupling El . 
Namely, the rotational momentum J even in the ground 
state (|i| = 5) mediates the coupling between levels and 
brings in the nonlinearity due to the reflection symmetry 
0, EH EI • Consequently, excited spectra espe- 
cially at big \j I are marked by a high density of avoided 
level crossings 0. The source of both MAC and local- 
ization is the interplay of the quantum coupling of levels 
and nonlinearity. The asymmetry of the interaction con- 
stants a 7^ (3 in E®(&i+&2) JT model bears an additional 
source of the quantum non-integrability. 



Recently, Yamasaki et al. [l5| first investigated the E(g>e 
JT model in terms of a search for possible quantum 
chaotic patterns. The semiclassical decoupling was per- 
formed according to the lines of the adiabatic approach. 
The principal attention was paid however to an explicit 
nonlinear term |lq | respec ting the trigonal bulk symme- 
try. In the previous work |l7j we investigated the quan- 
tum E®e JT model numerically and analytically in sev- 
eral limiting cases. We showed that an accurate account 
of the hidden nonlinearity can lead to nontrivial patterns 
similar to those produced in a system of two nonlincarly 
coupled oscillators. Our numerical analysis showed the 
presence of the chaotic motion domain at intermediate 
values of energy, which reflected in MAC in the quantum 
spectrum. Appropriate patterns also revealed from the 
statistical analysis of the nearest neighbor spacing dis- 
tribution (NNS) and the distribution of the "level curva- 
tures" (second derivatives of energies with respect to the 
coupling parameter a). 

In the present paper we investigate the generalized 
model assuming a ^ (3 [E<£>(&i + b 2 ) model]. The lo- 
cal spinless double degenerate electron level linearly cou- 
pled to two intramolecular vibron (phonon) modes is de- 
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FIG. 1: Complex behavior of energy levels of the nonsym- 
metric JT model as function of /3 - avoided crossings, crossings 
and clusterings of levels. 
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FIG. 2: (Color online) Examples of wave functions in the 
plane (Qi <X> Q2) for a = l,/3 = 0.9. Pronounced exotic states 
close to either displaced Fock states or "radial" states of the 
symmetric problem are shown. The extended (delocalized) 
state n = 190 is shown for comparison. 

scribed by the Hamiltonian 

H = {b\b t + b\b 2 + l)I + a{b\+b x )a z - (3(b\ + b 2 )(T x , (1) 

where a x , a z are 2x2 Pauli matrices, / is a unit ma- 
trix and the pseudospin notation refers to the two-level 
electron system. The operators 6j, b\ satisfy boson com- 
mutation rules = Sij. The interaction term oc a 
removes the electron degeneracy and the term ex (3 me- 
diates the phonon-assisted tunnelling. 

The Hamiltonian Q has SU(2) symmetry and com- 
mutes with the reflection (parity) operator 

R = R ph <Tx, R P h = exp(iirb\bi), (2) 

with RphQi = -QiRph, Qi = b\ + bi\ thus the eigen- 
states of the problem are chosen to have a definite parity 
p = ±1. The phase plane (a, (3) includes two limiting ef- 
fectively one-parameter models of higher symmetry: (i) 
the rotation symmetric limit a = /3, i.e. E®e Jahn- Teller 
case investigated previously an d (h) the one-phonon 
(Holstein) model of either (3 = or a = 0. In the latter 
cases the wave functions are coherent (displaced) Fock 
states exp(7(&j — bi))\n). 

For arbitrary (a, 0) the Hamiltonian can be exactly 
diagonalized in the electron subspace using the Fulton- 
Gouterman (FG) unitary operator ^lj 

u =W:X) (3) 

In the radial coordinates Q\ — > rcoscj) and Q 2 — > rsui(f) 




FIG. 3: (Color online) Effective potential built on coherent 
states of both phonons. Its minima visualize complex inter- 
play of effective quantum oscillators. 

the FG transformed Hamiltonian Q for p = +1 is writ- 
ten as [l3 

6 S mtl " = -Yri (4) + V - + 
V%ar (cos (f>I — sin (j)R p h^z) + v2(a — (3)r sin (f)R p h&z (4) 

[here the reflection (J2J acts as R p h(r, <j>)f(r, 4>) — f(r,TT — 
(/))]. Investigation of Etgie case a — (3 in terms of rota- 
tional quantum numbers [eigenvalues of conserved angu- 
lar momentum J — i(foi&J — b± b 2 ) — l/2g„1 has a long 
history, dating back to the paper in [18j- The spectrum 
separates into irreducible representations, each charac- 
terizing by the quantum number = i, |, The 

whole matrix is block-diagonal, and switching the term 
cx (a — (3) causes a complex interference of the levels in 
different blocks. The complicated structure of the en- 
ergy levels is exemplified in Fig. 1 where we show 40 
subsequent excited energy levels for a = 2 and varying 
f3. The level avoidings are accompanied by level degen- 
eracies (crossings). As in the symmetric case, there is 
a number of excited wave functions showing anomalous 
localization (Fig. 2) called "exotic states" pj- One can 
recognize two distinct types of the "exotic" states which 
are remnants of said limiting models with higher symme- 
try: in Fig. 2 the state n = 185 reminds us of the coherent 
state for (3 = while the states n = 69 and 186 with their 
markedly radial structure exemplify wave functions typ- 
ical for the case a — (3. 

Far from the rotation symmetry it is more convenient 
to perform the transformation © in the space Q\ x Q 2 : 

Hfg = b l b * + 1 + a ( b i +bl ^- P^ b l + b ^ R P h ■ ( 5 ) 

i=l,2 

The elimination of electron degrees of freedom reveals the 
nonlinearity hidden in the initial Hamiltonian Q [terms 
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with R p h in 01 or 0- Hamiltonian JSJ) differs from that 
of exciton (dimer) by the phonon-2 assistance in the tun- 
neling term + b 2 )R p h which accounts for Rabi os- 
cillations by the virtual emission and absorption of the 
phonon-1. These oscillations are essentially the origin of 
the nonlinearity of the reflection symmetric model and 
of its quantum nature. Namely, a consecutively classical 
version of the model would require to set R p h = ±1 drop- 
ping the nonlinear term. Thus the equivalence between 
the models and JHJ is lost and the classical analog of 
the model Q with (3 ^ is self-controversial. 

Averaging the diagonalized Hamiltonian over the 
trial wave function chosen as a combination of coherent 
states exp[i J2 Rph(jkPk + KkQk)]\0) maps it onto a 

fc=l,2 

Husimi form [n| for two oscillators nonlinearly coupled 
in variables 7,, TTj. The "effective potential" (Husimi rep- 
resentation of the Hamiltonian operator) built on the co- 
herent states in terms of the classical coordinates 71 , 72 
(Fig. 3) was used for a variational treatment of the 
ground state problem 0] . We showed that the parame- 
ter space (a, [3) is separated onto two regions: the region 
a > (3 of the dominating self-trapping or "heavy po- 
laron" with large 71 ~ —a and small 72 ~ /3exp(— 2^/\) 
and that (a < [3) of the dominating tunneling between 
electron levels or " light polaron" solution of small 71 and 
large 72. Quantum fluctuations cause mixing at the bor- 
der between them and thus broaden the transition re- 
gion. The potential in Fig. 3 visualizes the complex 
interplay of two oscillator potential wells resulting in the 
emerging of the third very narrow local minimum (at 
71 ~ and 72 > 0) responsible for the appearance of 
the light polaron (for more detailed discussion on this 
point see [13>LIJ])- This heuristic visualization can be a 
guide for understanding phenomena in the excited spec- 
trum as well. Excited states will follow the structure of 
either dominating self-trapping (a) or tunneling {(3) in- 
teractions. 

The level spacing distributions are considered as lead- 
ing characteristics to distinguish between quantum in- 
tegrability and chaos, the latter being described by the 
Wigner surmise (WD) P W (S) = (n/2)S exp(-7r5 2 /4), 
and the former by the Poisson statistics (PD) Pp(S) = 
cxp(— 5) [TCj, |2(J, |2l|] . Since long ago the Wigner surmise 
was conjectured as a limiting case for quantum chaotic 
behavior and supported from the point of view of ran- 
dom matrix theory (RMT; it is an exact conjecture for 
the 2x2 RMT version, and a rather close approximation 
to the exactly solvable case of random matrices of infinite 
dimensions 20] ) . The PD is associated with the superpo- 
sition of the multitude of uncorrelated levels: the falling 
exponential is a limiting case of a big number of inde- 
pendent level sequences, irres pec tively to the level distri- 
butions inside each sequence |22|. An adequate random- 
matrix model for our case of broken symmetry requires 
separately treating block and interblock elements |23ll24| 
yielding complicated statistical predictions as to the re- 
sulting superposition. Numerous interpolation formulas 
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FIG. 4: (Color online) Nearest-neighbor distributions of lev- 
els (unfolded and scaled to (S) = 1) for different values of 
a, (3. The curves on the histograms represent Poisson (long 
dashed), Wigner (full), semi- Poisson (stars) and yS (crosses) 
distributions. 



describing the intermediate situations between the com- 
plete integrability and chaos were considered [ToL l23| . A 
simple interpolation based on the information theory con- 
siderations was suggested [2j| in the form of the superpo- 
sition of WD and PD: P(S) ~ (7r/2)5exp(-^i5 - vS 2 ). 
The accommodation constants fj, and v had to be cho- 
sen to ensure the normalization with (S) = I conditions 
and yielding some a priori given variance cr 2 , thus giv- 
ing a one-parameter interpolation because the variance 
monotonically changes from (4/-7T — 1) (WD) to 1 (PD). 
Later on the theoretical background for this form of in- 
terpolation was enhanced on the base of the stochastic 
reformulation of the level statistics problem (2(|. The 
semi- Poisson distribution (sPD) P s p(S) = 45exp(— 25) 
with the dispersion exactly 0.5 pertains to this class. Re- 
cently it was introduced to mimic new seemingly univer- 
sal properties in certain classes of systems, in particular, 
being characteristics of the "critical quantum chaos" [27J , 
therefore it is worthy to probe it as third reference point. 

Figures 4 and 5 give the presentation of NNS statistics 
of JT system with varying a and (3. The level spacings 
were unfolded according to the common procedure in or- 
der to exclude the secular changes of level density and to 
ensure (5) = 1. Figure 4 gives examples of NNS distri- 
butions for sample values of a and j3. Figure 5 presents 
the standard deviations cr 2 = (5 2 ) — (5) 2 as functions 
of (a, (3). In the vicinity of either a = /9, or q > /3, 
a <C f3, the distribution functions P(S) follow rather 
Poisson statistics (Fig. 4, first row), but small deviations 
of parameters from these cases abruptly bring the dis- 
tributions away from it. Figure 5, however, reveals an 
interesting opposite universality for the parameters far 
from the one-parametric cases (that is supposedly in the 
most chaotic domain). The distribution functions seem 
to tend to a well-defined limit, but this limiting case does 
not resemble the Wigner surmise as one would expect. 
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FIG. 5: Standard deviations a 2 = (S 2 ) - 1 of the NNS distri- 
butions. Grid lines correspond to semi-Poisson (0.5) and \/~S 
(§) values. 



The standard deviation of all curves has lower boundary 
equal to 0.5 (in Fig. 5 it is markedly seen for the curves 
corresponding to a > 2 in the domains approximately 
0.4a < (3 < 2.5a but apart from (3 ~ a). Note also 
in passing a remarkable mirror symmetry of interchange 
a <-> f3: the corresponding level statistics are identical 
to a high degree of accuracy (the interchanged models 
however are not equivalent, and, for example, the wave- 
functions of the ground state for a > [3 and a < [3 differ 
essentially [II G3). 

Thus, the concurrence of two phonon modes of the 
Hamiltonian and the existence of two symmetry-changing 
limits bring into being a new limiting distribution which 
can be considered as a "most quantum chaotic" one for 
this system. We cannot conclude whether this hypothetic 
distribution is universal in the sense that it can encounter 
elsewhere. In the scope of the present paper we only try 
to mimic it discovering its possible universal properties. 
Natural suggestion is to compare it with the semi-Poisson 
distribution sharing the same cr 2 = 0.5 (Fig. 4, second 
row). Quantitatively, a coefficient of deviation from uni- 
versal limits of Pw{S) and Pp{S) can be introduced in 

the form V = J Q S °[P(S) - P w (S)]dS/ J So [P P (S) - 
Pw(S)]dS which gives the weight of the distribution left 
to the point So ~ 0.4729 of the intersection of WD and 
PD and ranges from (WD) to 1 (PD). For sPD rj sp ~ 
0.3858. The graphs of -q similar to that of Fig. 5 again 
show the existence of a lower boundary 77 ~ 4.7 > r) sp - 
Hence, the sPD turns out to be much a better fit than 
WD, especially at S > 1. A systematic shift of the mass 
of the distribution to the left with respect to sPD is 
markedly observed in Fig. 4. This shift is accounted for 
the nonlinear level repulsion of the actual distributions. 
At small S < 1 they scale as P ~ S s . The actual re- 
pulsion index (Brody parameter) S for "chaotic sets" of 
(a, (3) varies between ~ 0.3 and 0.5, meanwhile the class 
of WD and sPD assumes a linear repulsion P(S) ~ S. 
Seeking for universal properties of a limiting distribution 



we use the maximum value S = 0.5 and suggest another 
trial form P sq (S) = (3^/3/V^)y/Sexp(-3S/2). The co- 
efficient in the exponent is chosen to ensure (S) — 1 and 
to conform with the behavior of sPD at large S. Samples 
in Fig. 4 (second row) show that P SI? -distribution reason- 
ably fits actual distributions for small S in the chaotic do- 
main, although the latter have a tendency to shift slightly 
to the right. Therefore, the NNS distributions far from 
the E ® e and Holstein limits appear to be confined be- 
tween two suggested fitting formulas (note correspond- 
ing grid lines in Fig. 5). The standard x 2 reliability test 
however shows significant deviations between actual dis- 
tributions and both sPD and P sq indicating that neither 
of the reference distributions is a good fit in the strictly 
statistical sense. 

The sPD was recently suggested to describe a narrow 
intermediate region between insulating and conducting 
regimes exemplified by the Anderson localization model 
|27| , the mentioned opposite cases being described by cor- 
respondingly Poisson and Wigner statistics. At present 
a plausible analytical support (in the sense of RMT ap- 
proaches) for this new distribution and for its universal 
character is lacking. It was found numerically |28j that 
the width of this intermediate domain strongly depends 
on the length L of the system (its number of sites): for 
a system of infinite length one would get sPD in the 
narrow region of the Anderson parameter around W cr , 
meanwhile for short len gths the width of the intermedi- 
ate domain widens (in (23 the lengths of the order of 
~ 5 — 10 were checked). J aim- Teller systems from this 
point of view can be considered as systems with L = 2 
(two electronic levels regarded as pseudo-sites); thus it is 
to expect that the " transition domain" for such a system 
is rather large and plain in the parameter space, from 
whence there follows the applicability of sPD for almost 
all reasonable values of parameters a, (3. From the same 
point of view the systematic shift of the distribution to 
the left of sPD indicates that the generalized JT system is 
always closer to the " insulator" phase meaning the dom- 
inance of "heavy" polarons localized on one electronic 
level rather than the "light" ones. It is not surprising 
since heavy polarons are associated with the broad well 
of the effective potential (Fig. 3) which has a higher den- 
sity of states than the narrow one responsible for light 
polarons. On the other hand, the empirically s ugg ested 
second reference distribution is of Brody type [lfj with 
the scaling P ~ \fS at small S. The Brody parameter 
can be related to the fraction of chaotic motion ar- 
eas in the semiclassical picture. The present work gives 
merely a sketch of these possible relations, whose quan- 
titative examination would be a challenge for the future 
study. 
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